A method for the evaluation of the epsilon expansion of multi-loop masslessFeynman integrals is introduced. This method is based on the Gegenbauerpolynomial technique and the expansion of the Gamma function in terms ofharmonic sums. Algorithms for the evaluation of nested and harmonic sums areused to reduce the expressions to get analytical or numerical results for theexpansion coefficients. Methods to increase the precision of numerical resultsare discussed.
展开▼